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We present an elementary introduction to some of the quark models used to un¬ 
derstand the properties of light mesons and baryons. These lectures are intended 
for both theoretical and experimental graduate students beginning their study of 
the strong interaction. 


1 Introduction 

Modern nuclear and particle physics has a secure theoretical foundation based 
on the idea of local gauge invariance, known as the standard model. It has 
satisfied every experimental or theoretical challenge directed at it so far. In 
the electroweak sector it was spectacularly confirmed by the discovery of the 
heavy vector bosons (the and Z°) in the early 80’s. On the other hand, 
in the strong interaction sector, quantum chromodynamics (QCD) has pre¬ 
sented more problems. While it has passed every test in the regions where 
it can be solved (notably the QCD evolution of parton distributions in deep 
inelastic scattering), it has not been possible to solve it satisfactorily in the 
non-perturbative sector, especially for the structure and properties of the light 
hadrons, such as the nucleon. 

The fundamental particles upon which QCD is built are the quarks. In 
the early 1960’s there was a rapid proliferation in the number of “elementary” 
strongly interacting particles (hadrons) and it was suggested that these were 
more likely composite particles. The more fundamental particles from which 
they were built were initially called Aces or Quarks. By choosing them as a 
fundamental representation of the unitary symmetry group SU(3), one could 
immediately bring order to the chaos that had previously reigned in hadron 
spectroscopy. One unsatisfactory feature of the quark hypothesis was that no- 
one had observed a fractionally charged particle. They were sought in every 
conceivable place, from deep ocean sediments to rocks brought back from the 
moon. It well known that none were ever detected - at least as free particles. 
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On the other hand, there has been tremendous progress on the strong in¬ 
teractions, both theoretically and experimentally. There are firm indications 
that non-perturbative QCD will never allow us to see free quarks, that is, that 
they are forever confined to the interior of hadronic systems in a colour singlet 
state. Nevertheless, there is an overwhelming amount of evidence to support 
the idea that quarks, with precisely the expected electroweak properties, do 
exist in the interior of the observed hadrons. Indeed, one of the most exciting 
challenges facing nuclear physics is to investigate the role played by quarks 
in explaining the properties of nuclear matter, from the normal densities of 
observed nuclei through to the much higher densities at the centre of neu¬ 
tron stars and the possible deconfinement transition in relativistic heavy ion 
collisions. 

Initially, only three types of quark were required to understand the known 
hadrons, the up, down and strange quarks {u,d,s). However, the discovery 
of the J/ij) particle in 1975 increased this list to include the charm quark, c. 
Later the discoveries of the T led to the inclusion of the bottom/beauty quark, 
b, and finally the top/truth quark, t. 

We shall be concerned with just the three light quarks which are of most 
interest in nuclear physics. Moreover, because of the space limitations, it is 
not possible to discuss all of the models which have been invented to represent 
the non-perturbative regime of QCD. Over time these models have tended to 
become more and more complicated and, in some cases very abstract, so that 
the beauty and simplicity which led to the quark model in the first place can 
often be lost. In these lectures we shall concentrate on just a few, relatively 
simple models: 

1. Non-relativistic Quark Models 

• centre of mass can be done exactly, 

• chiral symmetry - problems, 

• relativity - problems. 

2. Nature of confinement 

• vacuum structure, 

• soliton. 

3. Relativistic, confining, potential model. 

4. MIT Bag 

• symmetries. 
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• conserved currents. 

5. Chiral Bags. 

Furthermore, our presentation will be at a very elementary level. It is intended 
primarily for beginning graduate students interested in the strong interaction 
- from both the experimental and theoretical perspectives. Although set at 
a relatively elementary level, we hope that the beginning student will find 
enough background, not often presented now, to be able to better appreciate 
the current literature together with a few insights that may help to make the 
subject come alive. 


2 Harmonic Oscillator Shell Model 


The first model that we shall look at is the harmonic oscillator shell model. 
This is quite a simple model, but it is nevertheless helpful in understanding 
many features of the hadron spectrum, the model is built on the valence quark 
picture in which a hadron consists of just three, confined, quarks. The concept 
of confinement is represented by linking the quarks together by virtual springs, 
so that a baryon, for example, is described mathematically by the following 
Hamiltonian 


^ = E 


^ 1 




( 1 ) 


The harmonic oscillator has the convenient property that the centre of mass 
motion can be exactly separated from the internal dynamics. In addition, the 
internal structure can be written as an effective two-body problem, with two 
quarks combined as a single subsystem. This involves the following definitions 


A = ^ (rj -I- f 2 - 2r3), 

V6 


( 2 ) 

( 3 ) 


and finally 

(4) 

where the coordinates can be visualised as in Fig. I. 

These definitions can be used to simplify the form of the Hamiltonian, Eq. 
(1), so it is easier to visualise the two body nature. 


H = Hcm 


2m + 2m + 2 ^ ^ 


( 5 ) 
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Figure 1: Picture of Harmonic Oscillator Shell Model. 


where Pp{p\) is the momentum conjugate to p (A) and i?cM is the kinetic 
energy of the centre of mass motion of the three-quark system as a whole. 
Clearly the three-body Hamiltonian has separated into two distinct oscillators 
describing the internal structure of the hadron. In fact, in this equal mass 
case, the oscillator frequencies of both two-body systems are the same 


UJp=LOx = 



and defining a constant a, as 


( 6 ) 


a = (7) 

we find that we can solve for the following wavefunctions. 


Quantum No.’s 

Wavefunction 

Energy 

o 

II 

V’oo = 

Eq = -1- ^hto 

o 71? 

II II 

i’ll = J^p3^/m(/5)exp"“"(^"+^") 

El — 

El = “1“ ^hio 


etc... 



A physical model is of limited use if it cannot be compared with real data. 
Within the standard model, the quark masses are currently free parameters. 
In the oscillator model these are often set to the values (constituent quark 
masses): 

mu ^rud^ 340MeV, . . 

ms ~ 500MeV. '' '' 
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Having quarks of unequal mass slightly complicates the kinematics. Suppose 
we have two quarks of equal mass and one different, then we can set: 


mi = m 2 = m ; m 3 = m', 


so that, through the dehnitions 


mp = m ; m\ = 


3mm' 
( 2 m + m') ’ 


(9) 


( 10 ) 


we can assign masses to the two body subsystems. Using these definitions in the 
appropriate Hamiltonian we hnd that in the centre of mass system {Hem = 0): 


—* 9 —* 2 

jj PP I PA 
2 mp 2m\ 


3k 


(P^ + A^) 


( 11 ) 


and we notice that cOp ^ to 

As a simple consequence of the fact that ujp i 0 \, as noted in the previ¬ 
ous paragraph, the excitation energies of the p and A degrees of freedom are 
different. Even this simple observation has profound phenomenological con¬ 
sequences. One of the mysteries of the baryon spectrum is that there are far 
fewer states observed than predicted by the quark model. A possible expla¬ 
nation of this, advanced by Isgur and Karl can be illustrated by the usual 
experimental method for exciting baryon resonances. For example, in the 
K + N^H^K + N reaction we hnd that the uj\ degree of freedom tends 
to be excited rather than cOp. (This is a direct consequence of the different 
energies required to excite the strange and non-strange degrees of freedom.) 
As a result one will tend to be experimentally blind to those strange baryon 
resonances where the excitation energy has gone into the LOp coordinate. In or¬ 
der to test this idea one needs new ways of exciting baryon resonances and this 
is precisely what is planned at new, high-duty factor machines like CEBAF. 


u • ^-T-► • d 


• s 

Figure 2: Idealised picture of, for instance, a A. 
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3 Spectroscopy 

The hadron wavefunction is dependent upon the intrinsic properties of the 
quarks: 

= dt(space, spin, flavour, colour). 

For baryons it must be totally symmetric in the first three of these, space, spin, 
and flavour, because the colour wave function must be totally anti-symmetric 
in order to produce a colour singlet object. Taking the A as an example, we 
know that the space, spin, and flavour components of the wave function are 
each symmetric (S), 

A: S,S,S 

However for the N each of the spin and flavour components is of mixed sym¬ 
metry (MS), with a symmetric spatial wave function: 

N : S, MS, MS. 

To first order the strong interaction we expect SU(6) spin-flavour symmetry to 
be respected for the light quarks, so given one hadronic state we can generate 
many others with raising and lowering operators. Below we present an example 
(and some exercises for the reader) of how operating on the wavefunction of 
one particle predicts the existence of another. (The solutions to some of the 
exercises are included in the appendix.) 


Example — Spin-Isospin 

Take for example the spin-isospin (SI) wavefunction of the 


SI 


T T T 

U[, Ufl 


We can calculate the wavefunction of the |A+, = -|-|) 

with the lowering operator /_(oc ti_ -I- T 2 - -I- T3-) with 


ri_ |ui) = |fii), 


and 


= 0 . 

Therefore, after this operation we have 




M T, M t, d T), 


( 12 ) 

by operating on it 

(13) 

(14) 

(15) 
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where |~) = sum over all different permutations. 

Similarly, and this is left as an excise for the reader, we can see the change 
in the spin of the particle by operating with the spin lowering operator S -, 


A+,J, 



oc S- 


A+, 



M T, M T, d i) + 

Ex#l 


U T, M i, fi t)) • 


(16) 


As a hint for this exercise, the first wavefunction on the right hand side of the 
equality in Eq. (17) has 3 terms, and the second has 6 terms. The spin-up 
proton is the orthogonal state with the same values of I 3 and Sz, and has 
wavefunction 


\P T)si - (2 


M T, u t, d i) - 


■u T, M i, d T 


Ex #2 (17) 


For Exercise 3 find the SI wavefunction of |E° and [At). (Hint: Start 
from |S*+, Jz = -|-|) = u t) M t, s T^). 

Exercise 4: Non-relativistically the magnetic moment of a particle of 
charge g, and mass m is 

,7=^?. (18) 

It is suggested that the reader shows that 


{p \p T) ^ _3 

( expt. 2.79 "j 

V “ -1.91 d 


(19) 


Energy Levels 

Calculations in this model lead to expectations of finding the L = \, negative 
parity, states at hoj and = 1 (Is) or L = 2 (Od) states at 2Tioo. Nature 
however presents us with two surprising results: 

• The Roper, Pn (1450), which would naively be a Is excitation, occurs 
below the lowest negative parity “nucleon excited states” - the D 13 and 
Sii at ~1550 MeV. 

• There exists 200MeV between the Is and Od “2/ia;” states N(1680) F 15 . 
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In fact the Roper is still a mystery. In the bag models it has been described as 
a “breathing mode” but it has also been described as the result of coupling 
to the inelastic two-pion channels 

R Ntttt 
N rj 


4 One Gluon Exchange 

We have so far considered only the simplest shell model picture of baryon 
structure, with the quarks moving in a mean confining field. However, one 
expects that there should be some residual interaction which, motivated by 
QCD, is usually taken to be the non-relativistic reduction of the One Gluon 
Exchange diagram shown in Fig. 3. This is proportional to the product of two 
quark-gluon vertices, ^ 1^2 = Ai • A 2 . Using the fact that the eigenvalue of 
the total colour wave function for a baryon or meson must be zero, one easily 
finds: 


Ai • A2 ) — - (Ai • A2 

Baryons 2 \ / Mesons 


Ex #5 
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Hint : = -T 


( 20 ) 

( 21 ) 


m 


m 





g 


gluon. 


Naturally this quantity is colour invariant, even though the A are the colour 
matrices. Defining the strong coupling constant in the usual way 


El 

47r ’ 


a. 


( 22 ) 



we find that the one gluon exchange potential is given by 


VoGE 


2 



1 


4^2j,3 


- - -^S{r) - \ Si2 - • ^2S{r) 

r 6m^ 

((cti + 2(72) • (rx Pi) - ((72 + 2(7i) • (r X ^ 2 )) • 


(23) 


The terms inside the square brackets can be thought of as the Coulomb, Darwin 
and tensor (S '12 = (3(7i • fd 2 • r — Si • S 2 )) terms respectively, while the final 
term {Si ■ S 2 ) is a hyperfine interaction. 


Example #2 

For the A, all qq pairs have S' = 1, which implies that 

{Si ■ Sj}^ = +1 (24) 

For the N, there exists an equal probability of finding S = 1 and S = 3 pairs, 
so we find that 

(<?i • (7 j)jy = ^ (+1 - 3) =-1 (25) 

Therefore, when the hyperfine interaction is treated as a perturbation, toa in¬ 
creases and ttin decreases, by equal amounts. This effect breaks the degeneracy 
of the naive model. 

One should note that on dimensional grounds 

(-ffhyp)(7V or A) (26) 

where the difference toa — toat determines the value of o^, which is typically 
of the order 0.6. 


Example :^3 The S(1180) A(1115) mass difference 


For the S we have 


(E| {Su ■ Sd |S) = -1-1, 
(S| Ss ■ {Su + Sd |S) = -4, 


and therefore 

, OsO' 

(^hyp/s OC ^2 

Similarly for the A we find 


dOsO^ 

fhms 


(E| {Su ■ Sd |S) = -3, 

(E| Ss ■ {Su + Sd |S) = 0, 


(27) 

(28) 


(29) 
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(i?hyp)ACX-^. (30) 

Hence we see that uiy: > toa- 

5 Hadronic Shell Model 

A slightly more sophisticated approach is the hadronic shell model, suggested 
by Isgur and Karl*. Defining the Hamiltonian as 

^ = ^ + + H^yP + Hs,. ( 31 ) 

^<j L J 

with q;^ ~ 0.6, we see that the l/r^ term is a colour Coulomb interaction 
while b represents a linear confining potential of strength ~ O.lSGeV^. The 
next step is to diagonalize H in, say, a 2Tilo space of harmonic oscillator wave- 
functions. This diagonalisation allows the calculation of baryon energy levels 
in the model, however it is usually done with the caveat that the overall scale 
must be adjusted for each major shell. The detailed wave functions given by 
the model allow one to investigate various electric and magnetic transition 
probabilities in detail. 



For example, for the neutron, [dd) pairs are always S' = 1, (which feel a 
repulsive hyperfine interaction) while the {ud) pairs are 3:1, (S = 0) : (S = 1) 
(attractive : repulsive). Therefore the effect of diagonalisation is to force the 
dd pairs further apart. 
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Open questions 


After studying these models one should be asking, among others, the following 
questions: 

• What is the origin of m? At the IGeV scale in QCD we have 

m„ ~ 6MeV, 

md ~ 9MeV, (34) 

rris ~ 150MeV. 

• What is the nature of confinement? For heavy quarks lattice QCD sug¬ 
gests that: 

E{R)r^-- + cr, (35) 

r 

captures the long and short distance pieces of the interaction. On the 
other hand, phenomenologically the spacing between levels in heavy 
quark systems {bb and cc) is more or less independent of the quark mass, 
which favours a fractional power dependence, like 

• Is it really a few-body problem? This question is directly related to the 
nature of the vacuum. 


6 Soliton Models 

One alternative to the models we have been discussing is the soliton model 
suggested by Lee®. Suppose that the vacuum solution of QCD is a colour 
dia-electric {k < 1). To quote T. D. Lee: 

“Quark confinement is a large scale phenomenon. Therefore, 
at least on phenomenological level, it should be understandable 
through a quasi-classical microscopic theory.” 

To implement this concept, suppose that there exists a medium with Kmed ^ 1- 
In such a medium we find that any charge will produce a hole (vacuum) with 
K = 1 (as illustrated in Fig. 4). 

It can be shown that the repulsion of s and the surface charges necessarily 
means that work has to be done to shrink the hole. The continuity of the 
displacement, D: 


Din(— Ain) — Dout(— ^medAout) 


r—R 


(36) 
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Figure 4: The dia-electric is anti-shielding. Here Kmed 1- 


means that the colour electric field is given by 


Bin 

Bout 




r < R, 
r > R. 


Thus we find that 


(37) 


{Electric field energy of cavity — Electric energy without medium effect} 

= / l-E ■ D — {Electric energy w/o medium effect} (38) 

a All Space ^ 

Prove, as an exercise, that the above expression has the form 
(^-1 _ 

- L^oo aS itmed ^ 0, Ex #8 (39) 

K 

Hence we can see that a perfect dia-electric is confining. 

Now we suppose that the dia-electric vacuum is a lower energy state, so 
that it costs energy to make this hole. The energy cost of making this hole is 
given by 

Ui^oie = BV + CS (40) 

where V is the volume of the hole and S is the surface energy. Therefore we 
find the equilibrium radius, iJequih when 


LL , 

i^hole ^electric) 

an 


= 0 . 


-R—^equil 


(41) 


Naturally we have that i?equii ^ oo when K^ed —*■ 0. Therefore there are 
only solutions of the form illustrated by Fig. 5, where the total charge inside 
the cavity is zero. This figure illustrates the close analogy between the super¬ 
conducting state, which excludes magnetic helds and the QCD vacuum which 
excludes the colour electric field, leading to confinement. 
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Figure 5: Solution form for the soliton model. Here the dashed lines can be 
visualised as the self-consistent “bag surface” 


Practical Implementation 

These ideas were implemented in the model known as the Friedberg-Lee Soliton 
6,7 

Cfl = qiip - TO - fa)q - ^kG^‘'-U{ a) (42) 

where ct is a new scalar field, and the term involving k is a dielectric function 
that we treat perturbatively for colourless states. 



Figure 6: Illustration of the scalar field and quark-scalar-field contributions 
to the total energy for two choices of the valence quark, scalar density. 


In Fig. 6 we see that the lowest energy state is when <7 = 0, which is the 
non-perturbative vacuum with no valence quarks. 

fqq)i (7 « 0 which means that when the scalar, valence quark density 
is high the perturbative vacuum is restored. 

fqq )2 CT « (Tvac which means that when the quark field is small the 
non-perturbative vacuum is restored. 

The results of this model are as follows 

• The quarks dig a self-consistent “hole” in the scalar held. 
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Figure 7: Formation of a non-topological soliton as a self-consistent combina¬ 
tion of localised valence quark density and a “hole” in the non-perturbative 
vacuum field configuration. 


• Asymptotically quarks have effective mass for fa^ac i+m) (where m is 
the perturbative QCD mass). 

• Quarks are only “confined” in the sense that if, for example, ~ (1 — 
-2—)" then K ^ 0 outside the soliton. Hence the total colour electric 

•^vac ' _^ 

energy is infinite, unless ^ = 0 — that is, it is a colourless state. 

This is shown by using the same argument as for a charge in a dia-electric 
cavity. 

7 Colour Dielectric Model 

Another interesting model is the Colour Dielectric ModeP’®d°.This model gives 
us an effective description of the long-distance effects in QCD. Once again, we 
can formulate a Lagrangian to allow the investigation of the theory: 


C = -u(x) + - K(x)Tr(G^j.G'"'"). (43) 

X 2 ^ 

Here we find that ctvX is a confining scalar field that mimics the non- 
perturbative effect of the gluons. Note that when x ^ 1 (inside the soliton) 
the second term goes to the usual mass term, rrnpip, while as x 0 the effective 
mass of the quark goes to infinity, which certainly implies confinement: 

TYl 

X ^ 0 ^-> oo. (44) 

X 

In the simplest case, the potential term, U{x), may be taken to be quadratic 
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in x: 


(45) 


u{x) 




X 


2 


8 MIT Bag Model 

The MIT bag model is of a similar vintage to the harmonic oscillator shell 
model discussed earlier. We shall take some care in explaining it because it is 
still proving to be extremely valuable in theoretical nuclear physics. In partic¬ 
ular, a large number of applications of the quark model of nucleon structure 
have been developed using the MIT bag. Two that we mention especially 
are the calculation of nucleon properties in medium (for instance in a neutron 
star), and the Quark Meson Coupling (QMC) modelAs a relativistic model 
which permits analytic solution, the MIT bag also allows us to examine the role 
of dynamical chiral symmetry breaking in hadron structure - a topic of great 
current interest (for example) under the heading chiral perturbation theory. 

A model developed by Bogolubov^^ in the late 1960’s was the basis for the 
development of the MIT bag model. The model was an attempt to phenomeno¬ 
logically describe confined, relativistic quarks in a finite region of space. We 
can view this model as, for example, an analytically solvable version of the 
Colour Dielectric model (where here x is 0 outside and I inside the bag). Bo- 
golubov considered the simplest case: a Dirac particle of mass m, moving freely 
within a spherical volume of radius R, in a scalar potential. 

Vs{r) =—9{R —r)m (46) 

Inside the potential, the effective mass of the quark is 0, and outside the 
potential (or bag) it is infinite, thus confining the quark inside the bag. 

Much of the structure of this model relies upon the time independent Dirac 
equation 

(47) 

where the Hamiltonian® is 


H = d ■ ■[> + P{Vs{r) + m). 


(48) 


The eigenstates of the Hamiltonian are classified by finding operators that 
commute with H. One of these operators is 


3 = 1 + 


<^(4x4) 

2 


“See Appendix 10 for the mathematical conventions used in this paper. 


(49) 
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with a necessarily of the form 


CT(4x4) 



(50) 


The other operator is the relativistic analog of an operator involving both spin 
and orbital angular momentum: 

k = a ■ I + (51) 

The relativistic analog is the obvious generalisation 

= (3{a-l + l). (52) 

So we can construct eigenstates ot I + s thus 

^ ^ = IXk) , (53) 

and hence define 

^|X^) = -«|XC)- (54) 

Using the above definitions, and the identities listed in the Appendix, it can 

be shown that H does indeed commute with both j and K, and that j and K 
commute with each other, that is 


[H, j] = 0 = [H, K] = [j, K]. Ex #9 (55) 


The easiest way to obtain the eigenvalues of K is to square it, and then operate 
on a wave function ijj. So we find: 


and hence 




( 3 ^ 


{a ■ -\- 2a ■ I + 1 


= 1^ + a - l + l 


= T 


1 

4’ 




(? + 

{j{j + 1 ) + \)'<P 

(j + 


(56) 


(57) 
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We see that the eigenvalues of K are k = ±(j + i), and thus 


(58) 


I for j = I — i 

—I — 1 for j = I + 2 


thereby dehning both I and j by the one quantum number k. 
Now, using 

^ , d 


we can show that 


Substituting this back into Eq. (47) and solving we hnd that f/' has a general 
solution of the form 

9{r)x'i, 


- i- X Ex #10 

r r 

(59) 

0 i 

id • f— —^—a • r (Bk — 1). 
or r 

(60) 


U/(r)xO„ 

which satisfies the following coupled, ordinary differential equations 


(61) 


K.g 


{E - Vs{r) - m)g = - 
{E + V,{r)+m)f= 1 + i I + 

For the most elementary case, k = —1 (sp, these equations simplify to 


^ + L 

dr r 
dr r 


Kf 

•) 

r 


dr"^ 


+ {E^ - (m + Vsy)u = 0, 


(62) 


(63) 


where the substitution u = rg has been made. Naturally there is a similar 
equation for I (= rf). It is left as an exercise for the reader to show that 
by following the method of conhnement, suggested by Bogolubov which 
involves dehning the scalar potential as 


T,. . f —m for r < R 
= | oforr>i?, 

it can be show that, requiring u and I continuous at r = R implies 


/ N \/(l — (E/mY) , , sinCER) 

cos{ER) + — -—-—^ sin(Ei?) = 


1 + (E/m) 


ER 


1 - 


E 


E + 

Ex #11 


(64) 


(65) 
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This gives us a boundary condition for this model. By taking the confining 
limit (m ^ oo) we get the eigenvalue condition, given by the relationship 
between the spherical Bessel functions 

jo{ER)=n{ER). (66) 


The energy levels may be parameterised by the definition 

En. = ^, (67) 

where n is the principle quantum number. We now find that the first two 
positive roots of the k = —1 states are: 


ni_i = 2.04 <—> Isi, 
^2-1 = 5.40 <—> 2si, 

and further solutions are easily computed. 

Thus Eq. (61) can be rewritten as 


= V'n.-l(r) 


= 


IG 


70 m 
rji m 


X-i- 


( 68 ) 

(69) 


(70) 




Figure 8: Pictorially we can see that u(r) = l(r) at the boundary (r = R) 


Confinement is achieved in the MIT bag model by requiring that there is 
no quark current flow through the surface of the bag. The (charge) density for 
the Dirac equation is given by 


f = P = V'7°V' 


OC 



(71) 
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so that the charge density is proportional to the baryon number density. We 
also note that the space components of the current are given by 


j{r) = ip^ip = (72) 

so it is found that 

= ijoji - ijoji = 0. (74) 

In particular, at r = R, we can see that r • j = 0, so there is no flow of 
quark current through the surface. 

As a side issue, the reader might like to note that for a typical hadronic 
radius, i? ~ 1 fm, the charge radius for a proton is 

(r2);f = 0.82fm, (75) 

and 

^ ~ 400MeV, (76) 

R 

which should be compared with the constituent quark mass. 


Magnetic moment 

We now consider the case of a charged quark in this potential and introduce a 
constant magnetic field B to the system. Naturally we have 

1 = Is X r. Ex #12 (77) 

The quark current has a magnetic moment jl that satisfies 


—fl ■ B = 



(78) 


By making a substitution for j (Eq. 72), and for A (Eq. 77) we have 


—jl • B 



dV ■ (^B X rj 

dVrx (^ip^d^jj) ■ B. 


(79) 
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Noting that both sides involve a dot product with B we find 


where 



/ 


dvfx ■4!'''dip 


: (algebra) 

— Mconf^; 


eR 

dconi = ^ 


2 rriconst 


412-3 

6(12-1) 

{0.83} 


(80) 


(81) 

(82) 


and as noted earlier, mconst = 12/i? is like a constituent quark mass. One 
should note that the previous relationship for i^(= —|) is preserved. 


Axial current 


We know that axial current is very important for weak interactions, and that 
it is calculated by 


where r is the 3 Pauli matrices for isospin and 


(83) 


(84) 


where both u and d are a four component spinor. 

In the non-relativistic limit Ai dominates. The gamma matrices give us 


7°775 = 


1 0 
0 -1 


0 a 

—a 0 


0 1 
1 0 


a 0 
0 a 


(85) 


and in the non-relativistic situation we have nearly, but not exactly equal 
l' 

, so 


to 


with equality when xp = 


A^ ~ 


- the index i is an index in isospin. 


(86) 
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Another exercise for the reader is to show, for the non-relativistic quark 
model, using the spin-isospin wave function obtained earlier, that: 


(pTIEL 


\p T) 


(Ptl^b T) 


5 

3’ 


Ex #13 


(87) 


where the index j labels the quarks, and the denominator involves simply the 
baryon spin-isospin operators. Thus, with the spin-flavor wave functions we 
have used and the nucleon axial charge operator, we find that 


9 ± ^ ( (pT|A|p t) \ ^ 5 


( 88 ) 


However, from the /3-decay of free neutrons we know that it is 1.26! 

One suggestion to remedy this problem with qa was that there might be a 
strong tensor force which would mix a sizeable d-state (L=2) component into 
the nucleon wavefunction, thus reducing qa- Exploration of the deformation of 
the A, through the E2/M1 ratio in the A ^ Ay transition has since suggested 
that this is unlikely to be the correct explanation - the N and A do not appear 
to have a significant, intrinsic deformation. 

On the other hand, for our case the quarks are not non-relativistic: 




{jo(^)} CT + (CT-r)CT(a-f)|ji 


T 

2' 


(89) 

(90) 


and 


J dH'i/'(r)775^V’(fO =A/'^ 





Ex #14 


(91) 


(For reference 




(numerical) ^ _»T 

= O.DOfT—. 
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This implies that 


„bag r 

^ = 0.65 X - 
5 v 3 

= 1.09, 


(92) 


a result which constitutes nearly a 50% improvement in the error and thus 
may be considered a major success of relativistic quantum mechanics. However 
this success is tempered with the knowledge that the quarks were taken to be 
massless. If we were to give the quarks a “current quark” mass, the lower 
component would decrease and hence gA increase, so that in the heavy quark 
limit gA increases to |. 

This model of Bogolubov produces many interesting results, but the value 
of the bag radius, R, is chosen in an ad hoc manner. What is needed is a 
self-consistent link between the size of the cavity, and what is in it, similar to 
the soliton models. This self-consistency was supplied by the MIT bag model 
in a fully covariant way. We specialise to the spherical, static cavity (which is 
the only case that is easily solved in 3-1-1 dimensions ^ ). 


Lagrangian formulation 

In the static cavity approximation the Lagrangian density describing the MIT 
bag is: 

Lca:v = [q{i^-m)q-B]eY-^qq5s, (93) 

where the last term results in the quarks being infinitely massive at the surface 
of the cavity. For the present we have not included the gluons for simplicity. 
(To include the gluons we must include a term —jTt{G • G), with 


J 1 inside the cavity 
[ 0 outside the cavity 

(94) 

and make the substitution d D.) We have defined: 


6 »v = 9{R-r), 

(95) 

h = 5{r - R), 

(96) 

and the Euler-Lagrange equation 


dC dC 

dq ^d{d^q) 

(97) 


^The deformed case has been looked at by Viollier, Kerman, and others 
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is satisfied. We demand that the action is stationary under the following 
transformations 


q + 6q, 
q^ q + 6q, 

R-^ R + 5R. (98) 

The final transformation results in 

0v + 6s{^R), 

5s ^ 5s — n- d5s, (99) 


where = (0, f) is defined to be an outward normal. 

Exercise 15 for the reader is to show that these transformations lead to 
the following equations 


i^q = mq r < R, 

Dirac Equation 
• nq = q r = R, 

linear boundary condition (l.b.c.) 

B = —in • d{qq), r = R, 

non-linear boundary condition (n.l.b.c.) 


^Ex#15 (100) 


Exercise 16 is to show that the linear boundary condition implies that 
at r = i?, 


= +qq 

= -qq 

= 0, Ex #16 (101) 


As qq vanishes at r = i?, it can be shown that the total energy of three 
identical quarks is given by 

EiR) = 3 ^ + ^R^B, ( 102 ) 

where the first term comes from the Dirac equation and the second term can 
be thought of as the volume energy. We also find the non-linear boundary 
condition is the same as the Lee soliton, that is 


dE{R) 

dR 


0 = AttBR? 


^ R = 



(103) 
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Thus we can choose B once and fix it, and then calculate the size of each of the 
hadrons since the radius is determined by what energy state the quarks are in. 
To get a instinctive feeling for this equation one can visualise it as requiring 
that the pressure from the vacuum is balanced by the energy of the quarks at 
the surface of the bag. 

Since this model, as described, is quite simplistic, possible improvements 
are easy to imagine and often easy to implement. One of the most important 
extensions was the inclusion of the zero point and centre of mass corrections 

E{R) ^ EiR) - (104) 

Also one gluon exchange may be included by solving Maxwell’s equation in a 
cavity, subject to the confining boundary conditions. This leads to a hyperfine 
interaction which has the same dependence on spin and colour as the hyperfine 
interaction discussed earlier and hence we get the required N — A, E — A mass 
splitting, etc. 

Given Ccav we can follow the common procedure used in QCD and look at 
the symmetries, and hence the conserved currents. Recall Noether’s theorem, 
that if is invariant (i.e. SC = 0) under the transformation 


where e is an infinitesimal constant, then 

dL 


f = 




h 


(105) 


(106) 


is conserved. Thus we have that if £ = £symm + £break(</>i)j then from the 
Euler-Lagrange equation: 


= 


^£break 

d(t>j 


fr 


We now will look at some of the symmetries in the model. 


(107) 


Example :^5 U(l) invariance 

It is easy to show that the Lagrangian density 

£cav = [iq^q - B]ev - ^qq6s, (108) 
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is invariant under the phase transformations 

q + ieq, 
q ^ q — ieq. 

Thus we find that the baryon number current (up to a factor 1/3) 

= {iqi^) {-iq)0w 

= qj^qOv, 

is conserved. 

Example Isospin Current 

If we define q = ^ ^ ^ ’ f that >Ccav is invariant under 

T • e 

<1 + 9 , 

T • S 

q ^ q-iq—^. (Ill) 

It is left as an exercise for the reader to check that the isospin current is 
conserved and is given by 

r = qi^^qOY. Ex #17 (112) 

Example #7 Axial Current 

For two-flavour, massless, QCD, C is invariant under 

T • e 

9 ^ 9-*^—759, 

T • S 

q^q-iq—^i5, (113) 

which are chiral transformations, and there exists a conserved axial current 

^^ = 97'"75^9- (114) 

However for C the surface term “—^qqSs” is not invariant under the chiral 
transformation, Eq. (113), and therefore 

= qj^-f5^q0v, (115) 


(109) 


( 110 ) 
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satisfies 


(116) 


T 

d^iA^ = -iq'y5-^qSs, 

and we see (as illustrated in Fig. 9) that the confining boundary violates chiral 
symmetry - i.e. mixes left- and right-handed quarks. This is a major problem 
since symmetries should always be a crucial guide in constructing models. 




Bag 

Wall 




Bag 

Wall 


Incident (Helicity +1) Reflected (Helicity -1) 

Figure 9: Violation of chiral symmetry at the bag surface 


But QCD already has a problem. If = 0 then 


/All Space 


dVdf,A>^ = 0 , 


(117) 


implies 


dVA^ >=- dVV • A = 0, 


(118) 


by the Gauss Theorem. By defining the axial charge as = J dVAP we find 
that it is a constant of the motion, that is 


[H,QP\=Q. (119) 

Thus for all positive parity eigenstates of the Hamiltonian there exists a de¬ 
generate, negative parity state, i.e. 

H |7V+) = m \N+) Qs \N+) = \N-) has 

H\N-)=m\N-). (120) 


This is clearly not seen in nature! 

The solution comes through the Goldstone Theorem. Either these de¬ 
generate, negative parity states exist OR |0) =/= 0. That is, as a result 
of spontaneous symmetry breaking, there exists massless, pseudoscalar Gold- 
stone bosons. As the first option is clearly incorrect, in an effective low energy 
description of QCD, we require a massless “pion” field, </>, in addition to the 
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quarks, so that >Ccav is invariant. Naturally, if we probe this at high momen¬ 
tum transfer the tt will show its internal qq structure. This is the basis of the 
Cloudy Bag Model discussed in Sec. 9, which has a Lagrangian of the form 

rcBM = [^#<7 -B]ev- ^qe^^-^'f^/fq6s + ^ {D^4>f , (121) 

where 

Df,4)= +f sin (122) 

= df,cP + 0{^^). (123) 

and £cbm is invariant under 

T • £ 

9 ^ 9 -*^— 759 , 

_T • e 

q^ q- iq^-j5, 

cj) ^ (l) +ef + f (e X (f)^ X $ l-ycot(y) . (124) 

Thus 

= 97^75^90v + + Yi9''(^sin(y) 

= Kn.rX + fd^4> + 0{^^). (125) 

9 The Cloudy Bag Model (CBM) 

The basic premise for this model is that we work perturbatively about MIT bag 
model solutions. This can be compared to the Skyrme models or the topological 
soliton models, where the pion field dominates the dynamics. The question 
remains as to whether they are equivalent. For QCD in 1-1-1 dimensions this 
equivalence can be demonstrated, but in 3-1-1 dimensions the belief is that the 
equivalence is broken. 

The CBM Lagrangian is defined to be 

/IcBM = {iq^q - mqq - B) 9v - ^qqSs - 

+ ( 126 ) 
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and here the terms involving m use the current quark mass to slightly (par¬ 
tially) break exact chiral symmetry. We also see that chiral symmetry requires 
(at least) linear coupling of pions to the bag through the r • </> term. That 
is, the CBM Lagrangian reduces to the MIT bag plus a free pion field with 
interactions between the two dictated by chiral symmetry. 


-^CBM — "^MIT Bag -^Free—-TT T -^int- 
Hence this implies that 

H = HmIT + Htt + iJint 

= ^ -b ^ Wfcaj.afc-b ^ ^P^aakV^°‘+ h.c. 

a k (y.^(3^k 


(127) 

(128) 
(129) 


where {k = k^i) and 


(3a _ ^ ^ 

“ 2 / 


dV exp S(r - R) (/3| q(f)^ 5 Tiq{r) |q;) , 


(130) 


with (/3| and |q;) being bag states. We stress that the pion-baryon couplings 
involve form factors which arise naturally as a consequence of the internal 
structure of the baryons and which can be calculated using the underlying 
quark model - be it a bag or something more sophisticated. 

Another exercise for the reader is, for a = |A^), and /? = \N'), using 


= K-i = 




2i?3(0- l)sm^H 


(131) 


to prove that 


^N'N 

k,i 


- (2u;fc)-^/" 


i n ji{kR) 

y (n-i) kR 


3 

SI {N'\ E 3'aif^a ‘ k\N)si 

a—1 


Ex #18 


Hence we see that using = 



|, and 


it can be shown 


u{k) = 


‘jjijkR) 

kR 



V 


N'N 

k.i 


i 

(2u;fc)i/2 



u{k) (A^'l Tid • k\N ), 


(132) 


(133) 


(134) 
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and comparing to the usual ttNN coupling, we find 

g-jrNN _ 5A 

2mN 2 / ■ 

With a little rearrangement, we hnd that 


fg-KNN = gArriN, 


( 135 ) 


(136) 


which is the Goldberger-Treiman relation. 

From this same C we can derive various coupling constants and form- 
factors, for instance, AA^tt, AAtt, EAtt, EEtt, SStt, etc. 


Example :^8 The physical/dressed nucleon 

We have included some examples that may assist in the understanding of this 


model here. The physical nucleon, 
tonian with eigenvalue toat: 


fV ), is the eigenstate of the CBM Hamil- 


H 


N) = Ton 


N 


(137) 


This should be compared with the eigenstate of the MIT bag piece of the 
Hamiltonian: 

Hmit \N) = I A), 

Emit |A) = |A), 

e.t.c. (138) 


where the one gluon exchange splittings, which may or may not be equal to 
the mass splittings of the observed baryons, is included in the terms. In 
first order perturbation theory: 

a\ = Vz 11 A) + V \an,) ] (139) 

1 rnj<s - Wk - rua I 

s^/^(l^v) + l:c.(^AA + y;c,(^—Ali>—\ ,140) 

( ^ -UJk ^ TON - TOA - Wfc J 

The convergence of this perturbative expansion can be proven to be extremely 
good as long as R is large (> 0.7fm) - i.e. the form-factor is soft - and as 
long as we restrict our model space to low-lying bag states (e.g. A, A[i?,...]) 
- similar to the 2fiio restriction in the non-relativistic oscillator model. 
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Example #9 Charge distribution of n 

Suppose we ignore spin-spin effects (small in the bag), then we find that the 
neutron bag has no charge distribution. On the other hand, the wave function 
of the physical neutron is: 



Here we see that |n) contributes nothing as does |n7r°). The coefficients | 
and i should be instantly recognizable as the Clebsch-Gordon coefficients for 
isospin coupling. Since there is only one term (|_p7r^)) that contributes (at 
this order) we can now see that the leading order effect in producing a non¬ 
zero charge distribution for the neutron is the (proton bag) -I- (pion cloud) 
component of the neutron wave function 

It is easy to see that, as the bag has a sharp surface, this simple model will 
produce a zero in the neutron charge density at the bag radius Detailed 
calculations of the neutron charge distribution are complicated by the need 
to project a state of good total momentum and to remove spurious centre 
of mass motion (which cannot be treated exactly as it was in the oscillator 
model). Nevertheless, after all the necessary corrections are made this feature 
persists to quite good accuracy. In Fig. 10, we see the results of a very recent 
calculation^®, which illustrate the importance of the current measurements of 
the neutron charge form factor at MIT Bates, TJNAF and Mainz. 


Example #10 Corrections to bag mass 

To lowest, non-trivial, order (second order in i/int) we have 


(0) I ^/nN-k 

rriN = my H- 


yy{k) 


nml Jq uJkimN - ujr - EN{k)) 


| 4 /^ a . kyyk) 

3 Timl Jo uJk{mN - ujk - Eyk)) 


(142) 


Since both these diagrams are attractive there is a sizable correction of 
-300 to -400 MeV to the total energy. 

For comparison we look at the A — A system. It can be show that Fig. 
11(a) and Fig. 12(b) are equal in magnitude and that Fig. 12(a) is small. So 
we find that 


f 2 r poo 

Re(mA) = m^^ + / dk 

37rm2 [Jo 


yvj{k) 

Wfe(TOA - Wfe - Eiq{k)) 
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Improved CBM neutron charge density 



Figure 10: Neutron charge density calculated in the CBM for several choices 
of bag radius. Note that the peak in the negative charge density always occurs 
very near the chosen value of R. 


+ 


N 


N 

(a) 


N 



N 


Figure 11: Pion self-energy corrections for the nucleon 


\ r jt , 

16 irml [Jo UkiruA - ujk - i?A(fc))J ’ 


(143) 


where V is the Principal value of the integral and in this case results in severe 
cancellations. The first integral term is represented by Fig. 12(a) and the 
second integral represents Fig. 12(b). In this case the attraction leads to a 
downward shift of 200 to 300 MeV! 

Hence we can see that 100 to 200 MeV of the N — A splitting comes from 
the pion self energy, thus we find that the value of as is even smaller. 
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+ 


A 


A 

(b) 


A 


Figure 12: Pion self-energy corrections for the A 


Observation 

If as ~ 0.3, as this analysis of the N-A splitting suggested, the spin-orbit 
problem in the baryon spectrum is essentially resolved. 

Observation ^^2 

Q) How can one ever seriously believe “shell model” type studies if effects of 
channel coupling are so large 

A) You can’t! Each resonance should be studied as a coupled channel scatter¬ 
ing problem... Until recently this has been seen as involving too much work! 
So far it has only been taken seriously in a limited number of cases - see 
also the discussion of Thomas and Milled®. 

10 Conclusion 

In this very limited space we have tried to summarise a great deal of physics 
that is not often taught. It is assumed that students of the strong interaction 
have somehow “picked up” the information. In this hnal, brief section we would 
like to make a few connections between the simple models and ideas presented 
here and some modern research topics. 

In the past few years it has been realized that a hyperhne interaction 
involving spin and isospin, rather than spin and colour (as for one gluon ex¬ 
change) offers some phenomenological improvements in fitting the baryon spec¬ 
trum The motivation for such an interaction is that there should be some 
residual short-distance remnant of the long range pion cloud required by chiral 
symmetry that we presented here using the CBM. Indeed, it has sometimes 
been argued that one should include only this interaction with no one gluon 
exchange term. 

There is currently great experimental activity and related theoretical in¬ 
terest in exploring the chiral and deconfinement transitions expected at high 
baryon density. One may ask whether these phase transitions are coincident or 
even whether they exist at all. Recent work from several different approaches 
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to QCE^®’^®, suggests that they may be coincident, occurring at 3-4 times 
normal nuclear matter density. Moreover, these calculations provide a rather 
unexpected motivation for a model of hadron structure like the bag. Below 
the critical density for the chiral/deconfining phase transition a region of space 
containing current quarks is unstable and thus at hnite, but low baryon density 
one expects to form colourless bubbles of chirally restored phase - just like the 
MIT bag. There is still an enormous amount of work to be done to move from 
this observation to a real understanding of the best way to model the nucleon, 
but this work has provided important insight. 

Chiral perturbation theory is currently enjoying tremendous popularity, 
with many fascinating examples of experimental interest being explored 
However, it has also been realized that the usual formulation using, for ex¬ 
ample, dimensional regularization presents some problems. In particular, the 
chiral corrections grow rapidly with the mass of the Goldstone boson, whereas 
one would physically expect it to be less important. The CBM has exactly 
the behaviour as chiral perturbation theory as m-n —>■ 0, but the corrections to 
hadron properties actually get smaller as the Goldstone boson mass increases 
because of the physically motivated momentum cut-off in the model - recall 
this cut-off is related to the internal quark structure of the hadron. There is a 
great deal of important work to be done to combine these two approaches to 
the same underlying physics - e.g. see Ref. 

With regard to the chiral structure of the nucleon we must mention the 
very exciting experimental results from Fermilab in which, for the hrst time, 
we have been given a direct view of the asymmetry between d and u quarks 
in the nucleon sea . It seems most likely that this asymmetry, which was 
anticipated on the basis of the CBM^^, is at least in part due to the pion cloud 
of the nucleon required by chiral symmetry 

We cannot finish a review such as this without mentioning the very fun¬ 
damental questions at the heart of nuclear theory, for which models such as 
those described here have an important role to play. In particular, we would 
very much like to understand the role that nucleon substructure plays in un¬ 
derstanding the properties of finite nuclei as well as matter at higher density. 
As just one example, the quark meson coupling model which is based on 
the MIT bag, offers a completely new saturation mechanism for nuclear mat¬ 
ter that seems very natural. It also offers a consistent theoretical framework 
for calculating the changes in hadron properties (masses, from factors and so 
on) in a nuclear medium. There is a great deal of experimental interest in 
these questions and we hope that these notes may be of some assistance in 
understanding what is being done and perhaps in contributing to it. 
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Appendix 

Dirac Matrices and Spinors 

Throughout this paper we have followed the conventions of Bjorken and Drell 
The gamma matrices satisfy the following condition : 

= V""- (144) 

Where the 7 matrices are defined by: 

7 ° = /3, (145) 

7 = (3d, (146) 


75 = 7^ = f7°7^7^7^ 

= -*7o7i7273, 

and their values in the Dirac representation are: 


7° = 


7 = 


75 = 


I 0 
0 -I 

0 a 

—a 0 

0 I 
I 0 


Obviously 7 ° is hermitian, 7 * is antihermitian, 7 I = I, and finally 

{75,7^} = 0 . 




(147) 

(148) 

(149) 

(150) 

(151) 


Selected Solutions for Exercises 

Exercise 1 


We find the spin-flavour wave function of the Jz = ^, in the following way. 
We act with the spin lowering operator on each quark as shown below. 


A+, = 


(^,)- 


u t M T d t 
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= ;^{kiwTfiT) + |uTMic^T) + |MTuTfii) + 
\uldU^) + \uUld^) + \uUUl) + 

|rf i u T T) + M T M i d T) + M T M T d i)} 

= 5 M t M T C? M T M i T)) • 

Exercise 2 

We want to find the wave function of 

\p t) = uud, 0 , 

which is totally symmetric in Spin-Flavour. We now find all totally symmetric 
combinations of these spin and flavour combinations. 

(u T M T d i)symm(gpi„_Flav) = (u T M T ^ i) + (« T ^ i « T) 

i M t U T) 

and 

(w T W i T)symm(Spin-Flav) = (« T U i d T) + (« i U T d T) 

We now antisymmeterise with respect to, say, spin (ti), and we get 

\P T) = (2 u'l u'l di^ - 

where 18 = 2^ x 3 -I- (—1)^ x 6 (the factors 3 and 6 indicate the number of 
terms in the kets respectively). Alternatively, we can simply write down the 
normalized combination of and u'\ m J, d orthogonal to the 

A+, Jz = ^ state. 

Exercise 3 

We shall just present the solutions here, as the reader should be able to obtain 
these results alone by now. 

|S0 T) = ^ |2(n t d T S i)s - (w T d i s T)s - (w i d T S T)s) 
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1 


|AT) = 

Exercise 4 

We know that 


^/T8 V 

1 r 


(2 


ut dt s 1} - 


ut d IsU- 


■u J, d t s T 


V6 V 


w T M i s T) - 


u i d T s T 


E ^ = do 1^(1). - ^(5,). - • 


So we find 


(d Tl E ^ b T) = -j^do (2(u T u T c? i)s - (m T M i 0^ T)sl ^(Su)z 


18 

-^(*Sd)z - |2(u T M T C? i)s - (m T U i d t)s) 


18 


4(('utut d i)s| Edi^ l(w T uf di)s) 

i 

-((uf ui d t)slEd»^ Kw T W i c^T)s)| • 


Therefore we find that 


(dtlE^b T) = § {4. (-). 1-3 +4. 3 

i ^ 

where the factors in each of the terms of the sum are the factor multiplying 
the ket, the factor multiplying the spin operator, the sum of the spins of 
the relevant particle, and the number of terms in the ket, respectively. This 
simplifies to 

dp = (p Tl Edi-^ \p t) = do 

i 

Similarly we find that the wavefunction of the spin-up neutron is given by 

\n t) (X I- \p t) 

Therefore we find that , and hence — = — 
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Exercise 5 


Beginning with the fact that 

^Ai + A 2 + As^ \N) = 0. 

We fully expand the term that is squared. We then use the fact that the 
particles are indistinguishable. Thus we note that, for instance Ai • A 2 = Ai • A 3 , 
so we find 

SA^ + 3Ai • A 2 = 0, 

or rearranging we get that 

Ai • A 2 = —A^. 

Now, we use the hint that is given 


(the sum of the index a is a sum from 1 to 8 - due to the fact that we are 
working in SU(3)). Substituting back be get 


Ai • A 5 


Baryons 


1 

2 

1 

2 

1 

2 



Mesons 




8 

3’ 


as required. 


Exercise 9 


The requirement is to show that [H,j] = 0 = [H, K] = [j, K], Since the work¬ 
ing is similar in each case we shall only present a result for the last commutator. 
We have that 


ilK] 


r ^ ^ T 

I + —,a ■ I 


where the proportionality comes from the fact that K is actually multiplied by 
7 o. Thus we find 


[j, K] oc 
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— o’jHi, Ij] + 2 [o’i, o'j]lj 
j^k 

= 0 , 


as required. 

Exercise 11 

To show that u and I are continuous at r = i? we need to find solutions for 
them sX r < R and r > R and then match at the boundary. For the case of 
r < R, Eq. (65) simplifies to 

u" + E^u = 0 

Solving this for u we find that u = ylsin(i?r). For r > R we have 

u" + {E^ - m^)u = 0 

with the general solution u = a +/3 exp“'^'®^“'"^'’. Requiring that 

that u is finite as the value of r increases (i.e. u < oo as r oo) we find this 
general answer simplifies to u = (3 exp“'^'®^“’”^'’. 

We now require that these solutions match at the boundary (r = R). Thus 
we have 


u{r = R) = Asm{ER) 

=/Jexp-'/®™^, 

so we see the solution has the form 

u = Asm{ER) exp-^'®'-™'(’’-«) 

Now we require that f(r) is continuous at r = R. In the simple case we are 
examining (k = — 1) we see / is given by 

f={E+{m + Vs))-^^ 
dr 

Again, using the substitution g = f we have that 

dg 1 / du u 
dr r \ dr r 
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Substituting our two solutions for u at the point r = i? we get 


— E cos{ER) — 
E 


sin{ER) \ 


R 




1 


E + m 


— \/— E"^ sm{ER) — 


sin{ER) \ 

R ) 


Rearranging, we find what we sought to find 

, , a/I — (E/mV . . sm(ER) 

cos{ER) H--- j— —— sin(£;i?) = 


1 + {E/m) 


ER 


1 - 


E 


E + m 
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